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Abstract. We present new data of the dwarf galaxy NGC 1569 at 450 pm, 850 pm and 1200/an taken with 
SCUBA at the JCMT and the bolometer array at the IRAM 30m telescope. After including data from IRAS 
at 12, 25, 60 and 100 pm, we have successfully fitted the dust grain population model of Desert et al. (1990) 
to the observed midinfrared-to-millimeter spectrum. The fit requires a combination of both large and very small 
grains exposed to a strong radiation held as well as an enhancement of the number of very small grains relative 
to the number of large grains. We interpret this as the consequence of large grain destruction due to shocks in 
the turbulent interstellar medium of NGC 1569. The contribution of polyaromatic hydrocarbons (PAH’s) is found 
to be negligible. Comparison of the dust emission maps with an HI map of similar resolution shows that both 
dust and molecular gas distributions peak close to the radio continuum maximum and at a minimum in the HI 
distribution. From a comparison of these three maps and assuming that the gas-to-dust mass ratio is the same 
everywhere, we estimate the ratio of molecular hydrogen column density to integrated CO intensity to be about 
25 - 30 times the local Galactic value. The gas-to-dust ratio is 1500 - 2900, about an order of magnitude higher 
than in the Solar Neighbourhood. 
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1. Introduction 


Dwarf galaxies characteristically have low metallicities 
and consequently low dust abundances. In dwarf galaxies, 
both dust properties and amounts differ from those in spi¬ 
ral galaxies, as for instance i mplied by a difference in IRAS 
colours (cf. Melisse & Israel 1994 ). However, actual dust 
abundances and dust composition are only poorly known 
for dwarf galaxies. The dust mass of a galaxy can be esti¬ 
mated reliably only if good (sub)mm measurements allow 
to constrain the amounts of relatively cold dust which may 
dominate the total dust mass with only a very limited con¬ 
tributions to the emission at infrared wavelengths. Such 
data are available for a limited number of galaxies, and 
are especially scarce for faint dwarf galaxies. 

NGC 1569 (Arp 210; VII Zw 16) is a nearby irregu¬ 
lar dwarf galaxy at a distance of 2.2 Mpc (Israel 1988, 
hereafter 188). It is a member of the low-galactic latitude 
IC 342/Maffei 1/Maffei 2/Dw 1 group, containing at least 
15 dwarf galaxies (Huchtmeier et al. 200C| ). NGC 1569 is 
presentl y in the aftermath of a massiv e burst of sta r for¬ 
mation (188; Israel & De Bruyn 1988| ; Waller 1991 ). Its 
present star formation rate, derived from the Ha luminos¬ 
ity, is 0.4 M 0 yr -1 (Waller 1991). In the recent past, this 


galaxy has experienced a starburst whic h sta rted about 
1 — 2 x 10 7 yr ago as estimated by Israel (188) and about 
1.5 x 10 7 yr ago determined by Vallenari & Homans ( 1996 ) 
from colour-magnitude diagrams. The end of the star- 
burst, about 5 Myr ago, can be well dated by a k ink in the 
synchrotron spectrum (Israel & De Bruyn 1988] ) and pho¬ 
tometric studies (Vallenari & Bomans 1996| , Greggio et 
al 1998 ). Although a small galaxy with neutral atomic 
hydro gen (HI) dimensions of 3 x 2 kpc (Israel fe Va n 
Driel 1990| ; Stil & Israel, in preparation, hereafter [5101 ), 
NGC 1569 contains two extremely compact luminous star 
cluste rs A a nd B (Abies 1971 ; Arp & Sandage 1985 ; Aloisi 
et al. 2001 and references therein) with bolometric lumi¬ 
nosities of order 10® Lq located in a deep HI minimum. A 
third such cluster is embedded in bright emission nebulos¬ 
ity (FT. Israel & W. Wamsteker, unpu blished ; Prada et 
al. 1994 ; ; listed as No. 10 by Hunter et al. 2001) coincident 
with the peak of the radio continuum distribution. 


The HI observations show the galaxy to be in solid- 
body rotation out to a radius of l 7 (0.64 kpc); however, 
much of the HI is in chaotic motion ( |SI0l| ) . To the east of 
NGC 1569, an apparently counterrotating HI cloud con¬ 
nected by an HI bridge to NGC 1569 is observed (Stil & 
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Israel 1998). The brightest HI with column densities N(HI) 
w 2 x 10 2U cm -2 , occurs in the three main peaks of an HI 
ridge in the northern half of the galaxy. About a third of 
the total mass of NGC 1569 resides in its neutral atomic 


hydrogen (188, SI01 ). Weak CO emission is found close 
to the radio continuum, between two HI maxima. It was 
mapped in the J=l-0 and J=2-l transitions by Greve et 
al. (1996). Comparison with the J=l-0 CO detection in a 
larger beam by Taylor et al. ( 1998| ) suggests that the maps 
by Greve et al. contain virtually all CO emission from the 
western half of NGC 1569. Aperture synthesis maps of the 
same two CO transitions with the IRAM interferometer 
show distinct CO c louds of sizes ranging from 40 to 100 pc 
(Taylor et al. 1999 ). As these maps recover only a quarter 
of the single-dish flux, weak and more diffusely distributed 
CO must be present. CO(3-2) observati ons s how a high 
<7=3-2/J=2-l ratio of 1.4 (Miihle et al. 2001 ) indicating 
a warm molecular gas phase. 

The far-infrared emission from NGC 1569, observed 
with IRAS, is remarkably strong for a dwarf galaxy; the 
continuum spectrum indicates that the dust in the galaxy 
is exposed to intense radiation helds (188). 

In this paper we extend the far-infrared spectrum of 
NGC 1569 to (sub)millimeter wavelengths by presenting 
new observations obtained with the IRAM and JCMT 
telescopes. These new observations allow us, for the first 
time, to determine the amount of dust in this galaxy and 
to study its properties in detail. 

2. Observations 

Observations at 850/zm and at 450 //in were made in 
October 1997 and again in March 2000 with the SCUBA 
camera at the JCMT on Mauna Kea (Hawaii). As the 
quality of the later observations greatly surpasses that of 
the earlier, we only discuss the March 2000 observations 
in this paper. SCUBA consists of two bolometer arrays of 
91 elements at 450 //m and 37 elements at 850 // m, bot h 
with a held of view of about 2.3' (Holland et al. 1999| ), 
thus covering the whole galaxy. Both wavelengths were 
observed simultaneously in jiggle mode. The total inte¬ 
gration times, half on source, half on sky, was 5 hrs. The 
images were chopped with a throw of 120 ” at a frequency 
of 7.8 Hz. 

We applied the standard reduction procedure: this in¬ 
cludes flat-fielding, removal of transient spikes, correction 
of atmospheric opacity (rg 5 o ~ 0 . 3 ;t 45 o ~ 1 . 1 ); pointing 
correction and sky-removal. The data were calibrated by 
observation of the standard source CRL 618 (Tg 50 = 11.2 
Jy, F 450 = 4.56 Jy) immediately before and after the 
NGC 1569 observations. From the calibrator images, we 
determined a beam-size (FWHM) of 15.3 x 15.6 " at 850 
//m. 

Observations at 1200/zm were made in December 1998 
at the IRAM 30m telescope on Pico Veleta (Spain), 
with the 19-channel bolometer array of the Max-Planck- 
Institut fiir Radioastronomie (MPIfR). Additional maps 
at 1200/zm were obtained in April 1999 and again in March 


2000 with the 37-channel bolometer of the MPIfR. The 
beam-size is 10.8". The observations were done on-the-fly, 
with a wobbler throw of 46". Opacities ranged from 0.1 
to 0.3. The data were reduced in a standard manner in¬ 
cluding baseline subtraction, spike removal and sky noise 
removal. The maps were calibrated by observation of the 
planets Mars and Uranus. 

3. Integrated flux-densities 

Maps of NGC 1569 at the three observed wavelengths are 
shown in Fig. 1 and 2. We have optimized the SCUBA 
dataset by excluding the two samples out of ten that were 
taken under poor atmospheric opacity conditions. The 
noise levels of the maps are 2.5 mJy/beam (1200 //m), 4 
mJy/beam (850 /zm) and 68 mJy/beam (450 /zm). 

We determined total galaxy flux-densities by integrat¬ 
ing the maps over increasingly larger areas until the cu¬ 
mulative flux-densities thus obtained converged to a fi¬ 
nal value. However, it was found that low-level emission 
from NGC 1569 extended over most of the SCUBA held 
of view rendering an accurate determination of the (sky) 
zerolevel questionable. This was not a problem in the 
larger IRAM 1200/zm held. For this reason, we htted the 
850/zm and 450/zm cumulative growthcurves to the scaled 
growthcurve at 1200 /zm, allowing for a zerolevel offset 
of the 850/zm and 450/zm growthcurves. The scaling of 
the 1200 /zm growthcurve was uniquely defined by the re¬ 
quirement that the difference between the growthcurves 
increases as the square of the radius, as expected if the 
SCUBA maps had a zerolevel offset. 

In this way, we obtained total flux-densities S 1200 = 
250±30 mjy, Sg 5 o = 410±45 mJy and S 450 = 1820±700 
mJy (Table 1). The errors quoted are formal errors based 
on map noise and including, at 450/zm and 850/zm, the 
error introduced by the zerolevel correction. The latter 
error is, however, small: From the goodness of the ht we 
estimate it to be about 4% at both 450/zm and 850/zm. 
growthcurves. The actual uncertainty is larger because 
of calibration uncertainties. In particular at 450 /zm, the 
telescope error pattern and relatively high atmospheric 
opacities conspire to produce a large uncertainty at this 
wavelength. We estimate the total error, including opac¬ 
ity correction, calibration, noise of the map and zerolevel 
correction to be about 50% at 450 /zm and 30% at both 
850 /zm and 1200/zm. 

Especially at the longer wavelengths of 850 /zm and 
1200 /zm, the broad-band flux-densities determined from 
the maps may contain non-negligible contributions by 
thermal free-free continuum emission and CO line emis¬ 
sion. The contribution of the CO line to 1200/zm broad¬ 
band emission can be estimated from the <7=2-1 CO mea¬ 


surements by Greve et al. (1996). Adding up all the emis¬ 


sion detected by them, we obtain 7co = 13.6 K km s -1 at 
230 GHz. Allowing for the presence of a similar amount 
of undetected CO emission in the remainder of NGC 1569 
and using a conversion factor between T m b and flux den¬ 
sity of 4.6 Jy/K and an adopted instrumental bandwidth 
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Fig. 1 . Map of NGC 1569 at 1200 pm, smoothed to a reso¬ 
lution of 13". Contour values are at levels -n, n, 2n, 3n 
where n = 5 mJy/beam. The lowest contour corresponds 
to 2a, where a is the noise level of the map. The posi¬ 
tions of the starclusters A and B are indicated, as well as 
the maximum of the CO distribution (Tayler et al. 1998) 
and the position of the two most promin ent H II regions, 
catalogued as number 2 and 7 by Waller ( 1991 ). 


of 50 GHz for the bolometer, we obtain a maximum of 
2.1 mJy contributed by the J= 2-1 CO line to the total 
emission in the 1200 /im band. The J= 3-2 CO emission, 
contributing to the 850/irn broadband emission, has been 
measured by Meier et al. ( [2001 ). At the CO peak in the 
map by Greve et al. ( 1996 ), they find the J=3-2/J=l-0 
ratio to be unity, characteristic of rather hot and dense 
molecular gas. From this result, we estimate a maximum 
contribution of 10 mJy to the emission in the 70 GHz 
wide SCUBA band at 850 /*m. We thus find that, at ei¬ 
ther wavelength, contributions by CO to the total emis¬ 
sion of NGC 1569 are unimportant, reflecting the general 
weakness of CO emission from dwarf galaxies. 

Consideration of the thermal free-free continuum con¬ 
tribution does not warrant a similar conclusion. Its inten¬ 
sity has been estimat ed fro m radio maps at 1.5, 5, 8.4 and 
15 GHz by Wilding ( 199C ). The derived total thermal ra¬ 
dio continuum flux-density at 1 GHz is 100 mJy, in very 
good agreement with the value of 97 mJy estimated by 
Israel & De Bruyn ( 1988 ) from reddening-corrected Ha 
measurements. Scaling this emission to higher frequen¬ 
cies (shorter wavelengths) as z/ -0 ' 1 , we derive considerable 
thermal free-free contributions to the total emission of 58, 
56 and 52 mJy at 1200, 850 and 450 pm (23%, 13% and 
3%) respectively. In Table 1, we have summarized these 
results, together with further photometric data from the 
literature. 


4. Nature of the dust emission 
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Fig. 2. Maps of NGC 1569 at 850 and 450 pm. Contour 
values are at levels -n, n, 2n, 3n ...., where n = 8 mJy/beam 
(850 pm) and n = 115 mJy/beam (450pm). The lowest 
contours correspond to 2a, where a is the noise level of 
the map. 


4.1. Observed emission spectrum 


The dust emission and extinction in our Galaxy can be 
modelled as the combined effect of three components 
(e.g. Desert et al. 1990 hereafter DBP90; Siebenmorgen 
& Krugel 1992 ). These components are: ( i ) large grains 
obeying a power-law size distribution (cf. Mathis et al. 
1977), {it) polyaromatic hydrocarbons (PAH’s) responsi¬ 
ble for mid-infrared continuum and spectral feature emis¬ 
sion, {Hi) very small grains (VSG’s) with sizes of a few 
nm. The large grains are in equilibrium with the radiation 
field and their emission F\ can be described by a modi¬ 
fied black-body curve following from the Kirchhoff-law: 
F\ = B(A,T)k\, where B(A,T) is the Planck function, A 
the wavelength, T the dust temperature and k\ the ex¬ 
tinction coefficient. In contrast, the VSG’s are not at all 
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Table 1 . NGC 1569 emission from UV to millimetre 


Table 2. Flux densities of the various dust components 



UV and optical flux density 

A 

(A) 

flux density 

10 -14 erg s 1 cm -1 A -1 

Ref. 

1500 

152 ± 48 

(1) 

1800 

208 ± 25 

(1) 

2200 

146 ± 25 

(1) 

2500 

100 ± 18 

(1) 

3650 

108 ± 10 

(2) 

4400 

101 ± 9 

(2) 

5500 

70 ± 6 

(2) 


A 

H 

Integrated flux density from 12pm to 1200/im 
Observed Thermal+CO Dust Ref. 

(Jy) (Jy) (Jy) 

12 

0.59 

— 

0.59± 0.06 

(3) 

25 

5.95 

— 

5.95 ± 0.6 

(3) 

60 

44.6 

— 

44.6 ± 4.5 

(3) 

100 

52.2 

— 

52.2 ± 5.2 

(3) 

155 

36 

— 

36 ± 11 

(4) 

450 

1.82 

0.052 

1.77 ± 0.88 

(5) 

850 

0.41 

0.066 

0.34 ± 0.10 

(5) 

1200 

0.25 

0.060 

0.19 ± 0.06 

(5) 


(1) 188, corrected for Galactic foreground extinction as de¬ 
scribed therein. 


(2) De Vaucouleurs et al. (1991), correc ted for Galactic fore¬ 
ground extinction based on E(B — V) (188) and a standard 
solar neighbourhood reddening law as in 188. 

(3) IRAS point source catalogue, flux densities are colour- 
corrected as descr ibed t herein. 

(4) Hunter et al. (1989) 

(5) This work 


in equilibrium. Their heat capacity is so limited that they 
can be heated to very high temperatures (T ss 1000 K) 
by the absorption of a single photon. However, as they 
cool down rapidly, an ensemble of small grains cannot be 
described by a single temperature but rather by a rela¬ 
tively broad temperature distribution. Consequently, its 
integrated emission spectrum is significantly broader than 
that of large grains with a limited temperature range. 
Because the small-grain temperature is on average higher, 
this spectrum peaks in the mid-infrared. The wavelength 
dependence of k\ is not well known, however, models of 
interstellar dust that are both astronomically realistic in 
composition and able to reproduce the observed dust ex¬ 
tinction and emission curves, predict (3 cs 2 for large 
grains (e.g. Draine & Lee 1984, Ossenkopf & Henning 


1994, Kriigel & Siebenmorgen 1994). Observations of ac¬ 


tively star-forming galaxies, for which the dust emission 
beyond about 60 pin can be well described by only one 
temperature component, confirm a value of [3 close to 2 
(e.g. Chini et al. 1992, Kriigel et al 1998). A mo re general 
result has been obtained by Dunne et al. ( 2001 ) who have 
derived on a statistical basis (3 cs 2 for a sample of nearby 
galaxies with 450 and 850 pm data. For small amorphous 
grains, a lower val ue of (3 ~ 1 can perhaps be e xpect ed 
(Seki & Yamamoto 1989, Tielens & Allamandola 1987). 


Model A: Cold dust Model B: Very small grains 


A 

(/nn) 

Shot 

(Jy) 

*5*warm 

(Jy) 

^cold 

(Jy) 

SVsG 

(Jy) 

^large grains 

(Jy) 

12 

0.68 

— 

— 

0.44 

- 

25 

5.7 

0.18 

— 

6.07 

0.16 

60 

1.9 

41.2 

— 

26.2 

18.1 

100 

0.46 

52.5 

— 

25.3 

26.6 

155 

0.11 

26.6 

0.02 

17.2 

14.1 

450 

— 

1.1 

0.70 

1.9 

0.54 

850 

— 

0.12 

0.26 

0.39 

0.05 

1200 

— 

0.03 

0.11 

0.16 

0.01 


The integrated spectrum from 12pm to 1200pm of 
NGC 1569 provides important clues for the origin of the 
emission and for the relative importance of each of the 
three components described above. Flux-densities at wave¬ 
lengths A < 100pm define a slope much less steep than 
expected for a modified black-body curve. In particular 
flux-densities at 12pm and 25pm are much higher than 
expected for such a curve. Moreover, the steep rise in in¬ 
tensity from 12 to 25 pm precludes any significant contri¬ 
bution by emission from PAH’s, as the broadband spec¬ 
trum of these decreases with wavelength. At the long- 
wavelength side, the corrected flux-densities between 450 
and 1200 pm are nominally proportional to about A“ 2 ' 5 , 
i.e. proportional to a modified blackbody emission with 
a wavelength dependence (3 < 1, more characteristic for 
small than for large dust grains, or typical for a superpo¬ 
sition of various large grain temperature components. 

In order to further analyze the dust population of 
NGC 1569, represented by the observed emission, we have 
modelled the spectrum for a few representative cases. We 
present here two particular cases: (i) emission dominated 
by large particles at different temperatures and (ii) a 
three-component fit in which the relative fractions of dif¬ 
ferent types of dust particles are varied. 

4.2. An abundance of cold dust? 

In principle, it is well possible to fit the emission spectrum 
over the entire observed range with large grains at various 
temperatures. As we have only seven observed points, we 
must limit ourselves to a model allowing for three differ¬ 
ent temperature components. Obviously, even better fits 
can be obtained by allowing for more temperature com¬ 
ponents, but there is no physical relevance to such an im¬ 
provement. In Fig. ^jwe show the result, adopting (3 = 2. 
Table 2 presents the flux densities of the different compo¬ 
nents. The data can be fit very satisfactorily, which is no 
surprise given the large number of free parameters (almost 
equal to the number of independent data points). As Fig. || 
shows, each component is essentially determined by two 
different measured points. However, in order to reproduce 
the specific features at the short and the long wavelength 
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Fig. 3. The mid-infrared to millimeter spectrum of 
NGC 1569, fitted with a three-temperature large-grain 
model. Crosses mark flux-densities from Table 1. The 
three dust components have temperatures of 105 K 
(dashed line), 34.5 K (dotted line) and 7 K (dashed dotted 
line). 


sides referred to in the previous section, the inclusion of 
both rather hot and rather cold dust components is un¬ 
avoidable. The temperature of the coolest component re¬ 
quired to fit the relatively high (sub)millimeter emission 
is very low, of the order of only 7 K. 

Moreover, because of the low dust emissivities at this 
temperature, a very large fraction of all dust would have 
to be so cold in order to explain the emission at 850 /tin 
and 1200/tm. From the large grain component fits to the 
NGC 1569 spectrum we may estimate the dust mass from 
each components 100/rm emission: 


M d 


D 2 S 


100 (im 


B ioo/i m (^O^lOO/xm 


(1) 


Here, B WOll m(T) is the Planck function at 100/irn. We 
have assumed a dust extinction coefficient (per dust mass) 
of kioo/xm = 63 cm 2 g” 1 . Extrapolating this value with 
k\ oc A -2 and adopting a gas-to-dust mass ratio of 150 
we obtain fci 2 oo^m = 0.0025 cm 2 g” 1 (per gas mass) in 
good agreement with other studies (Kriigel et al. 1990, 
Mezger et al. 1990, Draine & Lee |1984 ). The results for 
the three components shown in Fig. |3| are a small mass of 
hot dust, Md,h = 12.4 Mg , a significant mass of warm 
dust, Md, w = 3.1 x 10 4 M 0 and a very large mass for the 
very cold dust: Md, c = 1.14 x 10 6 M 0 . The total dust 
mass Md,tot = 1.17 x 10 6 M 0 would thus be completely 
(i.e. for 97%) dominated by the cold dust component. 

How large is the uncertainty in this mass? The dust 
mass depends on the adopted dust temperatures and 
wavelength dependence of the extinction coefficient, (3. 


The highest temperature of the cold dust that is just 
marginally consistent with the errors of the 450 - 1200 
/mi data is 11 K. This temperature is mainly determined 
by the slope of the spectrum between 450 - 1200 /mi. With 
this, we derive a total dust mass of 5 x 10 5 M 0 , about a 
factor of 2 lower than derived for a dust temperature of 7 
K. 

The actual dust mass depends strongly on the choice 
of (3. For instance, if (3 = 1, we could fit the dust emission 
longwards of 60 /m i with only one dust component at a 
temperature of about 45 K - no cold dust being required 
at all. The total dust mass would then be much lower, 
only about 10 4 M 0 . However, we do not believe (3 = 1 
to be a plausible choice, because nearby galaxies provide 
compelling evidence (see Sect. 4.1) that (3 ~ 2 for large 
grains. Thus, if we had f3 ~ 1 in NGC 1569, it would 
indicate that its dust properties are very different from 
other nearby galaxies, a possibility no different from the 
one discussed in Sect. 4.3. 

In order to explore the influence of a slightly differ¬ 
ent /3, which might still be consistent with the observa¬ 
tions presented in Sect. 4.1, we have fitted the spectrum 
of NGC 1569 with f3 = 1.8. The results are very similar 
to (3 = 2. The hot dust component is unchanged and the 
warm and cold component require temperatures of 37 and 
7 K, respectively. The total dust mass is slightly reduced 
to 8.7 x 10 5 M 0 . 

In conclusion, if (3 cs 2 as in all or most other nearby 
galaxies, then the dust mass necessary to explain the spec¬ 
trum of NGC 1569 within a multi-temperature model is 
at least 5 x 10 5 M 0 and, if we take the best-fit to the data, 
even higher, 1.2 x 10 6 M 0 . 


4.2.1. Where can the cold dust hide? 


Such a large amount of cold dust is a very unlikely state of 
affairs in a low-metallicity galaxy dominated by intense ra¬ 
diation fields. For instance, from the UV measurements in 
|I88| we estimate a mean radiation field at 1000A within 
NGC 1569 (assuming NGC 1569 to be spherical with a 
radius of 1 arcmin = 640 pc) of about 1 x 10” 16 erg 
sec” 1 cm” 2 Hz” 1 . As the Solar Neighbourhood is char¬ 
acterised by a value of ab out 1 x 10” 18 erg sec” 1 cm” 2 
Hz” 1 (Mezger et al. 1982), the NGC 1569 radiation held 
held is two orders of magnitude higher in good agreement 
with the strong thermal free-free and Ha emission from 
this post-starburst galaxy. Such high radiation held inten¬ 
sities suggest that dust temperatures should be a factor of 
100 1 / 5 = 2.5 (for (3 = 1) to 100 1 / 6 = 2.1 (for (3 = 2) higher 
than for the Solar Neighbourhood. The latter has to be 
about 10 K for silicates and 18 K for graphite (Mathis et 
al. 1983) which would imply a minimum dust temperature 
of 20 K in NGC 1569. The same conclusion is drawn from 
a consideration of nearby spiral galaxies. Notwithstanding 
their higher metallicity (more shielding) and lower radi¬ 
ation held intensities, the coolest component contribut¬ 
ing to their far-infrared/submillimetre spectra has a tern- 
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perat ure of typically 15 K - 2 0 K 
et al. 


199S, Israel et al. 1999 


NGC 3627: Sievers et al. |1994 ). The coldest dust compo- 


e.g. NGC 891: Alton 
M51: Guelin et al. |1995|, 


nent in a galaxy such as NGC 1569 with very little shield¬ 
ing against very strong radiation fields should clearly be 
warmer than this. 

We can estimate the global dust amount from an en¬ 
ergy budget consideration by comparing the absorbed to 
the bolometric radiation. Integrating the UV flux densi¬ 
ties from 1500 - 3650 A, the optical flux densities from 
3650 - 5500 A and the infrared-to-millimetre flux den¬ 
sities from 12 - 1200 g m we obtain T\jv = 2-9 x 10 —9 , 
-Font = 1-7 x 10 -9 and A dus t = 3.7 x 10~ 


9 —1 

a erg sec 


cm 


We have to make a distinction between the dust emis¬ 
sion from HII regions and the diffuse dust emission. The 
dust in HII regions is mainly heated by ionizing photons 
that are almost completely absorbed by the dust locally, 
practically independent of the dust amount. The diffuse 
dust, on the other hand, is mainly heated by nonionizing 
photons, and the amount of radiation absorbed is directly 
related to the dust opacity. Thus, in order to calculate 
the (diffuse) dust opacity, we have to subtract the dust 
heating by ionizing photon. We estimate the ionizing UV 
radiation from the Ha emission. Waller (1991) finds a to¬ 
tal Ha luminosity of Ah« = 4.5 x 10 4 ° erg sec -1 . This 
value is corrected for the Galactic foreground extinction 
of E(B — V) = 0.56 (188); we neglect any extinction inter¬ 
nal in the HII regions. The total dust emission originat¬ 
ing from ionizing radiation can be estimated as ~ 


M dust = 8.8 x 10 4 x t b X 


kpc 2 


( 2 ) 


lower because the extinction efficiency per mass is higher 
for VSG’s (DBP90). 

The above considerations do not exclude the presence 
of very dense regions of cold dust barely contributing to 
the extinction or overall emission. But where could such a 
dust component be? It would have to be in regions of high 
dust opacity in order to be shielded from the interstellar 
radiation field (ISRF) or be situated far away from stars. 
The latter is not the case: NGC 1569 is a small galaxy 
with an intense star-formation over the whole area where 
dust emission is observed. In fact, the maximum of the 
dust emission is very close to the large star clusters A 
and B (Fig. 1). Is the self-shielding of the dust enough to 
maintain the inner parts cold? The opacity, r, along the 
line of sight can be calculated as: 


S„ 


SIB V {T) 


(3) 


27.12 x L Ha = 1.2 x 10 42 erg sec -1 (Xu & Buat pV95| ). The 
corresponding flux is F^ t = L^ a /(AnD 2 ) = 2.0 x 10 -9 
erg sec -1 cm -2 . 

The fraction of radiation absorbed by the diffuse dust 
is (Fd us t - F™ n st )/(F VV + A op t + Adust - F ™ t ) = 0.27. We 
estimate the diffuse dust opacity using a simplified radia¬ 
tion transfer model with a slab geometry, assuming that 
dust and stars are homogeneously mixed. The dust mass 
derived in this geometry is higher than if the dust were in 
a foreground layer. We use an approximate formula (Xu & 
de Zotti 1989 ) for the dust absorption probability, further 
simplified by applying for the UV and optical radiation 
the extinction properties a t 2000 and 4300 A, respectively 
(see Lisenfeld, Volk & Xu 1996 ). With this we derive an 
opacity in the blue of tb = 0.16. Assuming that the dust 
has the same extinction properties as dust in our Galaxy 
(DBP90), we can estimate the dust mass (in Mg ) as 


where S u is the flux density coming from the solid angle 
9. For a dust temperature of 7 K (11 K) we derive from 
the peak flux density of 25 mJy/beam ri 2 oo/jm = 1-5 x 
10 -3 (5.8 x 10 -4 ). With rv/ri 2 oo/jm = 4 x 10 4 (Kramer 
et al. 1998 ) this gives an optical opacity of Ty = 58 (23). 
Maximum shielding could be achieved if this dust were in 
one big cloud illuminated from outside only. But also in 
this case, a rather thick outside layer of warm dust exposed 
to the ISRF will be present. In order to attenuate the 
ISRF of NGC 1569 by a factor of 100, which would make 
it comparable to the Galactic ISRF, a layer of Ty = 4.5 is 
needed. Thus, we can estimate the cold dust mass fraction 
by considering a sphere with diameter Ty = 58 (23) where 
an outer layer of Ty = 4.5 is made of warm dust and the 
inner part consists of cold dust. The mass ratio of cold 
dust to total dust is 60 % (T = 7 K), respectively 20 
% (T = 11 K). This is not enough to explain the dust 
emission spectrum for which 97 % (respectively 94 % for 
T = 11 K) of the dust mass must be cold. 

An alternative way to shield the dust could be molec¬ 
ular gas which is most likely abundant in the regions of 
maximum dust emission (see Sect. 5). However, molecular 
hydrogen is dissociated only by photons of energies above 
11 eV, corresponding to radiation shortwards of about 
1000 A so that the overall radiation field of NGC 1569 
does not get significantly damped. 

A final argument against a large amount of cold dust, is 
that if it were present, the resulting dust mass is of order 
one per cent of the gas mass in the same area, M ras ss 
9.2 x 10 7 M 0 (HJ 


SI01; see also Sect. 4.3.2). This metal- 


where A is the surface area of the galaxy. Assuming a 
circular area with radius 80" we obtain Md us t = 2.3 x 10 4 
Mg . The uncertainty in this estimate is, even taking into 
account unknown parameters like the geometry of dust 
and stellar distribution, certainly less than a factor of 2 
so that the dust mass derived is more than an order of 
magnitude lower than the cold dust necessary to explain 
the dust emission spectrum. If the dust contained more 
VSG’s than in our Galaxy, the derived dust mass would be 


poor galaxy would thus have an usually high dust-to-gas 
ratio. 

For these reasons, we reject any explanation of the far- 
infrared/submillimetre spectrum of NGC 1569 involving 
the presence of large amounts of cold dust. 

4.3. An abundance of very small grains? 

As an alternative, we have fitted the dust emission spec¬ 
trum of NGC 1569 with the model by DBP90. In this 
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Fig. 4. The mid-infrared to millimetre spectrum of 
NGC 1569, fitted with the dust model of DBP90| assuming 
dust to be exposed to a radiation field equal to the Solar 
Neighbourhood held scaled up by a factor 60. Emission 
by PAH’s is negligible and the number of VSG’s is in¬ 
creased by a factor of seven over the local values in DBP90. 
Crosses mark flux-densities from Table 1. 


model, the dust is described by the three components men¬ 
tioned above (PAH’s, VSG’s and large grains) exposed to 
a radiation held which can be varied in intensity. The 
maximum and minimum size, as well as the exponent of 
the size distribution of the various components can be var¬ 
ied. Here, we use those values for which DBP90 achieved 
the best-ht for the solar neighbourhood. The big grain 
size ranges from 15 and 110 nm and they are distributed 
in a power-law with exponent 2.9. The VSGs have sizes 
between 1.2 and 15 nm with a power-law exponent of 2.6. 
The corresponding values for PAHs are 0.4 to 1.2 nm and 
an exponent of 3. The exponent of the wavelength de¬ 
pendence of the extinction coefficient is (3 = 2 for large 
grains and f3 = 1 for VSGs. We have tested different radi¬ 
ation fields and allowed for different relative contributions 
of the three components. We obtain a good fit for dust 
exposed to a radiation field similar in spectral shape to 
the local Solar Neighbourhood field, but with an intensity 
sixty times higher, consistent with the radiation field of 
NGC 1569 estimated in Sect. 4.2. 


The resulting fit is shown in Fig. |j and the flux den¬ 
sities of the various components are listed in Tab. 2. Not 
unexpectedly, the abundance of PAH’s is found to be neg¬ 
ligible because of the stee p mid- infrared rise in the re¬ 
quired spectrum. Madden ( 2000 ) reaches a similar con¬ 
clusion from observations of the MIR/FIR spectrum of 
NGC 1569. However, the abundance of VSG’s with re¬ 
spect to the large grains must be increased by a factors of 
about seven over the Galactic abundance assumed in the 


DBP90 model. The large VSG contribution is required to 
successfully match both the mid-infrared and the submil- 
limetre/millimetre ends of the spectrum. The contribution 
of large grains is well-constrained by the requirement to 
fit the far-infrared data points. 

The conclusion that the abundance of VSG’s is en¬ 
hanced with respect to the solar neighbourhood is very 
robust and does not depend on the details of the mod¬ 
elling, as we confirmed by trying to fit the data with other 
radiation fields given in DBP90 and changing the size dis¬ 
tribution of the VSG’s. Acceptable - although not equally 
good - fits could be achieved for the radiation field around 
an 05 star and for an enhanced UV radiation held. In 
those case the abundance by VSG’s had to be increased 
by a factor of 2 - 3. 

Thus, there are two principal differences between our 
model fit of NGC 1569 and the model fits presented by 
DBP90 for the local Galactic environment. These are: (i) 
the almost total absence of PAH emission and (ii) a strong 
enhancement of the VSG contribution relative to that of 
the large grains by a factor of 2 - 7. Both factors point 
to a different evolution of dust grains under the different 
environmental conditions prevalent in NGC 1569. 


4.3.1. Dust grain processing 


The first result is not entirely unique to NGC 1569. Both 
Magellanic Clouds also suffer from significant PAH deple¬ 


tion but not as stro ngly as NGC 1569 (Schwering |1988 ; 
Sauvage et al. 1990| ). However, it is relevant to note that, 
although the metallicity of NGC 1569 (12 + log {O/H) = 


8.19, Kobulnicky & Skillman 1997 ) is in between those of 
theLMC and the SMC, the intensity of its radiation held is 
5-10 times higher. PAH depletion, in fact, appears to be 
a general characteristic of dwarf galaxies. In the compila¬ 
tion by Melisse & Israel (1994), NGC 1569 clearly has the 
highest /25//12 ratio (about 10), but it is followed by two 
more galaxies (NGC 3738 and IC 4662) with ratios of 7 
and 8 respectively. The mean f '25 //12 ratio of the remain¬ 
ing ten Im galaxies is 4.5 ± 0.5, whereas 15 larger galaxies 
classihed Sm have a mean ratio 2.9 ± 0.4. Large late-type 
spiral galaxies with less intense mean radiation helds and 
higher metallicities mostly have ratios 1.5 - 2. With the 
/25 //12 ratio as indicator, irregular dwarf galaxies clearly 
have dust contents signihcantly depleted in PAH’s as com¬ 
pared to spiral galaxies. Most likely, PAH depletion is a 
function of the incident UV radiation, i.e. of the global 
UV radiation intensit y and the amount of shielding avail¬ 
able (Puget & Leger 1989 ). In this respect, NGC 1569 
provides one of the more extreme environments by being 
metal-poor and presently in a (post) starburst phase. 

The magnitude of the increase in the VSG/large grain 
ratio with respect to that of the Solar Neighbourhood, 
suggests the amount of VSGs increase at the expense of 
large grains by extensive processing of the latter. Grain- 
grain collisions in slow shocks can efficiently destroy large 
grains by fragmenting and shattering and convert them 
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into small grains producing an excess of small particles as 
we believe to see in NGC 1569 (Borkowski & Dwek 


1995 


Jones et al. 1996 ). NGC 1569 is a post-starburst galaxy. A 
large fraction of its interstellar medium is in chaotic mo¬ 
tion with a rather high mean velocity dispersion of about 
20 km s -1 , and over a dozen radio supernova remnants can 
be identified in this small galaxy. Thus, shocks are neces¬ 
sarily present in its interstellar medium. In addition, the 
low metallicity of NGC 1569 may also hamper formation 
of large grains by inefficient accretion to begin with. 


4.3.2. Dust mass and gas-to-dust ratio 


The total dust mass derived from the adapted DBP90 


5. Distribution of dust, atomic and molecular gas 

The (sub)millimetre maps not only closely resemble one 
another, but also maps of centimetre-wavelength radio 
continuum emission (Israel & de Bruyn 1988, Wilding 


1990) and Ha (Waller 1991) images. In all these maps 


model is much less than that found for the model based 
on cold large grains discussed in the preceding section. 
The mass of the large grains is derived from eq. (jlj) with 
'S'ioo fim and T being the flux density at 100 /tin of large 
grains and their temperature. The mass of the VSG com¬ 
ponent is found from the relative mass contributions of 
VSG’s and large grains given by DBP90| (their Table 2), 
scaled up by the factor 7 that we have found for NGC 1569. 
Notwithstanding the substantial number increase of small 
grains with respect to large grains, the total dust mass 
Md.tot = 3.2 x 10 4 M 0 is still dominated by the large 
grains, with M & t lg = 2.1 x 10 4 M 0 (66%). The large- 
grain temperature required by the model fit ( Tlg = 32 
K) is close to the temperature (Tiras ~ 35 K) that, ig¬ 
noring everything else, we would obtain directly from the 
IRAS / 60//100 flux-density ratio. Partly because of this, 
the dust mass derived here is somewhat fortuitously not 
very different from the mass 2.8 x 10 4 M 0 we would naively 
obtain from the IRAS 60/nn and 100/iin data only. This 
situation is very different indeed from the one pertinent 
to spiral galaxies where IRAS-derived dust masses com¬ 
monly underestimate the total dust mass by an order of 
magnitude (Devereux & Young 199C ). The uncertainty in 
this mass derivation can be estimated by considering the 
masses derived for the full range of fits consistent with 
the errors. In no case a lower dust mass was derived. The 
highest dust mass derived (in a cases where a lower VSG 
abundance was needed) was 6.3 x 10 4 M 0 , i.e. not even 
a factor of 2 higher. The range in dust mass derived from 
the DBP90 model is thus 3.2 — 6.3 x 10 4 M 0 . 


emission peaks in the western part of the optical image 
of the galaxy, just west of ‘super’ star cluster A (see Fig. 
1). The 850//m/T200/un continuum peak is close to the 
discrete CO clouds mapped by Greve et al. ( 1996| ) and 
Taylor et al. ( 1 999| ) , but offset from these by about 10" 
to the east. In all maps major emission curves eastwards, 
north of star clusters A and B, following a string of opti¬ 
cally prominent HII regions (cf. Waller 1991). The overall 
extent of the (sub)millimetre continuum emission is very 
similar to the optical extent, and that at radio continuum 
wavelengths. 

We might suspect that contamination of the 
(sub)millimetre continuum maps by free-free emission 
from ionized gas causes their resemblance to radio con¬ 
tinuum and Ha maps. Indeed, in Sect. 3. we have shown 
that our maps at 1200/im and at 850/im contain non- 
negligible contributions of free-free emission. However, the 
fraction of global emission involved (25% and 13% re¬ 
spectively) is insufficient to explain the resemblance es¬ 
pecially at 850/tm. Nevertheless, the question thus arises 
whether the free-free emission emission is proportional to 
the dust emission over the entire maps. If this is not the 
case, specific structural details in our maps could indeed 
be due to local concentrations of free-free emission. In 
order to address these concerns, we have compared the 


We have summed HI maps of NGC 1569 (SI01) over a 
circular area with radius of 80" comprising all the flux in 
the 1200/tm map. We find an HI mass of Mm = 5.1 x 10' 
M 0 . For mass contributions Mm. = 0.35Mhi (Israel pV97| ) 
and Mm = 0.25 M gas , and a dust mass M ^ = 3.2 — 6.3 x 
10 4 M 0 , we find an overall gas-to-dust ratio M gas /Md ~ 
1500 — 2900 which is significantly higher than the Solar 
Neighbourhood ratio. The corresponding ratio Mj-n/Afd ~ 
1600 is in g ood a greement with the prediction by Lisenfeld 
& Ferrara ( 1998| ) that this ratio should exceed the Solar 
Neighbourhood value typically by an order of magnitude 
for metallicities 12+log(0/H) = 8.70 and 8.19 in the Milky 
Way and NGC 1569 respectively. 


Ha map by Waller (1991; his Fig. 3a) to our maps. The 
Ha contour values multiplied by factors of 5.1 x 1CP 14 
and 4.9 x 10” 14 should correspond to the free-free emis¬ 
sion contribution across NGC 1569 in our 1200/im and 
850/im maps. With these conversion factors, the Ha map 
shows that the free-free emission from ionized gas is on 
the whole reasonably well scaled with the thermal emis¬ 
sion from dust. A notable exception is the main emission 
peak at a(2000) = 4 ft 30 m 47 s , <5(2000) = +64°50'58", co¬ 
incident with the most prominent HII region complex of 
NGC 1569 (HII 2 in Fig. 1). The ionized gas should con¬ 
tribute about 12 mjy/beam to the peak flux densities in 
our maps. This corresponds to disproportionate fractions 
of about 50% and 20% at 1200/tm and 850/im respectively, 
about twice the global average. Below emission levels of 15 
- 20 mjy/beam (1200/im) and 40 - 50 mjy/beam (850/im) 
ionized gas and dust emission are roughly proportional, 
equal to the global average. 

We have convolved to the 16" JCMT 850/tm resolution 
not only the IRAM 1200/im map, but also the WSRT HI 
column density map ( |3I01 ) (Fig. |[). The (sub)millimetre 
emission is mostly just south of the HI ridge in NGC 1569. 
The main peak, marking the location of molecular clouds 
and bright HII regions, is located on a saddle point in 
the HI distribution, i.e. close to a local minimum in the 
HI ridge, between the two elongated major HI maxima. 
We may use the maps for a direct comparison of the 
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Fig. 5. Intensity at 850/im (contours) over HI column 
density (greyscales). Contour levels are —8, 8, 16, 24, ... 
mJy/beam. Grayscales increase linearly from 1 x 10 21 to 
7 x 10 21 cm -2 . The hatched circle indicates the size of the 
16” beam (FWHM) applicable to both maps. 




Fig. 6. Correlation between HI column density and in¬ 
tensity at 850/im (left) and 1200/tm (right). Small and 
medium-sized black dots correspond to 5". The medium¬ 
sized dots mark positions where the 21-cm radio contin¬ 
uum intensity is higher than 20 mJy/beam. The large dots 
are binned averages of the remaining points (small dots) 
and the straight line is a fit to these average points. 


(sub)millimetre emission intensity with the column den¬ 
sity of neutral hydrogen (HI). Fig. 0 displays the corre¬ 
lation between the latter and the 850/im and 1200/im in¬ 
tensities. The overwhelming majority of points (small dots 
in Fig. 0) follows a linear relationship between broadband 
intensity and HI column density. At these intensity lev¬ 
els free-free emission is a more or less constant fraction 
of dust emission, and at these wavelengths most of the 
continuum emission is contributed by very small grains. 
Thus, the linear relation reflects a fairly constant VSG- 
to-gas ratio in NGC 1569. The relatively small fraction 


of points that deviate strongly from the general correla¬ 
tion is contributed entirely by the regions of most intense 
emission in the SCUBA/IRAM maps, notably the peak 
just mentioned. 

In order to determine the general correlation, we 
first applied a cutoff for all points with more than 20 
mJy/beam in the 21-cm radio continuum map; this ef¬ 
fectively excludes all strongly deviating points. We then 
averaged the dust emission in 5 x 10 20 cm -2 intervals of HI 
column density. The resulting points are marked by large 
dots in Fig. 0. A fit to the points thus found, corrected 
for the free-free contribution, yields the relations /1200 = 
6.2 x 10~ 22 N m - 0.51 and / 850 = 2.0 x 10" 21 Ahi - 1.71, 
(Nhi in cm' 2 , I in mJy per 16 arcsec beam). Assuming a 
constant gas-to-dust mass ratio, and neglecting local tem¬ 
perature variations, the deviating points then seem to im¬ 
ply the presence of a large amount of undetected hydrogen, 
which should be in molecular form. This consideration 
allows us to estimate the molecular mass independently 
from uncertain CO-to-H 2 conversion factors (Israel 1997). 

The largest difference is found for the peak in the 
maps, with a mean intensity 46 mJy per 16 arcsec beam 
at 850/un and 22 mJy per beam at 1200/tm. Corrected 
for a mean free-free-contribution of 10 mJy/beam (some¬ 
what less than the free-free peak), these intensities imply 
total hydrogen column densities .ZVh = 2 x 10 22 cm -2 . The 
difference IVh — -/Vhi is the column density of hydrogen 
atoms bound in molecules. Thus, Fig. 0 suggests a mean 
molecular hydrogen column density A7 h 2 = 8 • 10 21 cm~ 2 . 

Division by the mean integrated CO intensity over this 
area, Iqo ~ 1-5 K km s -1 (Greve et al. 1996| ) yields 
an estimate for the actual CO-to-H 2 conversion factor 
A 1569 ~ 5.3 x 10 21 cm -2 (K km s -1 )” 1 . This result, about 
25 - 30 times the commonly adopted Galact ic val ue, is 
very close to that obtained by Greve et al. ( 199C| ), but 
rejected by them; it is a factor of thre e belo w the high 
but uncertain value estimated by Israel (1997) and about 
a factor of three higher than the estimate by Taylor et al. 
( 1999| ) which was derived from virial considerations and 
therefore should be considered a lower limit under the con¬ 
ditions prevalent in NGC 1569 (Israel 1997; 200C). 

It is interesting to note that the deviating points in 
Fig. implying the need of molecular gas, only occur for 
IVhi > 3 x 10 -21 cm -2 . This could be due to a threshold 
in A^hi necessary for the formation of molecular gas. Such 
a threshold exists in our Galaxy around Nhi ~ 5 x 1CP 20 
cm -2 (Federman et al. 1970) and reflects the balance be¬ 
tween molecular gas formation and destruction by pho¬ 
todissociation. A minimum gas column density is neces¬ 
sary in order to protect the molecular gas from the ra¬ 
diation field through self-shielding as well as shielding by 
dust grains. It can be expected that in NGC 1569 the 
threshold TVhi is higher, in agreement with our finding, 
because the low dust content decreases both the shield¬ 
ing and the molecular gas formation rate and because the 
high radiation field increase the molecular gas destruction 
rate. 
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As can be seen in Fig. ^], adding these molecular hy¬ 
drogen column densities to those of HI results in a total 
filling-in of the HI ridge minimum. It appears that the to¬ 
tal hydrogen distribution in NGC 1569 remain ridge-like, 
but with a smooth increase to a pronounced maximum 
coincident with the present continuum peak. 


6. Conclusions 


1. We present new maps of the dwarf galaxy NGC 1569 
at 450, 850 and 1200 /irn taken with SCUBA at the 
JCMT and the MPIfR bolometer array at the IRAM 
30m telescope. Integrated flux-densities at these wave¬ 
lengths may be compared to those at at 12, 25, 60 and 
100 p,m obtained earlier with IRAS. 

2. The steep rise in intensity from 12 to 25 fin l excludes a 
significant contribution from PAH’s, as the broadband 
spectrum of these decreases in wavelength. 

3. The (sub)millimetre flux densities are high compared 
to the flux densities at shorter wavelengths and the 
(sub)millimetre spectrum has a relatively shallow slope 
(~ A -2 ' 5 ). Such a spectral shape can be explained by 
the presence of a significant amount of cold dust. Fits 
to the observed spectrum by three-temperature dust 
models, however, require most dust to be at tempera¬ 
ture of only about 7 K. 

4. We do not favour this explanation. The intense radi¬ 
ation field and low metallicity of NGC 1569, imply¬ 
ing poor shielding, render it very unlikely that large 
amounts of dust at such low temperatures may exist 
in NGC 1569. We show that the high dust opacities 
necessary to shield a large fraction of the dust from this 
radiation field are not present. Even if it could, the re¬ 
sulting gas-to-dust mass ratio would have to be about 
100 , again an unusually low value for a low metallicity 
galaxy such as NGC 1569 (12 + log (O/H) = 8.19). 

5. Alternatively, the (sub)millimetre flux densities may 
be dominated by emission of VSG’s at various non¬ 
equilibrium temperatures. The combined emission is 
characterized by a wavelength dependence mimicking 
an extinction coefficient k\ oc A -1 . Using the dust 


model of PBP90| , we achieve good fits for dust exposed 
to radiation fields similar in spectral shape to the Solar 
Neighbourhood field but with sixty times higher inten¬ 
sity, as is appropriate for NGC 1569 (188). 

6 . The fits require VSG to large grain abundances to be 
enhanced by a factor of 7 as compared to those in Solar 
Neighbourhood. The precise value of the enhancement 
factor is slightly model dependent; use of different radi¬ 
ation field yielding less good, but still acceptable fit re¬ 
sults, yields slightly different factors. A robust conclu¬ 
sion is that the VSG enhancement factor in NGC 1569 
is 2 - 7. 

7. Although in these models much of the emission origi¬ 
nates from very small grains, virtually all of the mass 
still resides in the large grains. The gas-to-dust mass 
ratio is 1500 - 2900, about an order of magnitude 
higher than in the solar neighbourhood. 


8 . Both the dust and molecular gas distributions peak at 
a local minimum in the HI ridge. If gas-to-dust mass 
ratio are constant over NGC 1569, the lack of HI at 
the very peak of the dust emission indicates the pres¬ 
ence of a significant column of molecular gas. From 
the inferred molecular gas column density and the ob¬ 
served CO emission, we estimate the hydrogen column 
density to integrated CO intensity conversion factor 
X ~ 5 x 10 21 H 2 mol cm -2 (K km s -1 ) -1 , or about 
25 - 30 times the local Galactic value. 
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